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Abstract 

We propose to extract the charge information of high energy muons in very inchned 
extensive air showers by analyzing their relative lateral positions in the shower 
transverse plane. We calculate the muon lateral deviation under the geomagnetic 
field and compare it to dispersive deviations from other causes. By our criterion 
of resolvability, positive and negative muons with energies above 10^ GeV will be 
clearly separated into two lobes if the shower zenith angle is larger than 70° . Thus 
we suggest a possible approach to measure the //"'"/^^ ratio for high energy muons. 
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1 Introduction 



The study of cosmic rays with primary energies above 10^ GeV are typically 
based on the measurements of extensive air showers (EAS) that they initiate 
in the atmosphere. The ground detector array records the secondary particles 
produced in shower cascades, including photons, electrons (positrons), muons, 
and some hadrons. Then their arrival times and density profiles are used to 
infer the primary energy and composition of the incident cosmic ray particle 
yy, usually through comparison with simulated results. EAS events initiated 
by primary particles with energies above 10^°'^ GeV, the so called GZK cutoff 
energy [2], have already been reported Questions about the composition 
of such ultra high energy (UHE) cosmic ray particles are still open to investi- 
gations 
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Photons, electrons and positrons are the most numerous secondary particles 
in an EAS event. However, for very inclined showers which will be concerned 
in this article, these electromagnetic components would travel a long slant 
distance and are almost completely absorbed before they reach the ground [5J. 
On the other hand, muons are decay products of charged mesons in shower 
hadronic cascades. Most high energy muons survive their propagation through 
the slant atmospheric depth, during which they lose typically a few tens of 
GeV's energy |6]. These high energy muons carry important information about 
the nature of the primary cosmic ray hadron, which will be extracted from 
their energy spectrum and lateral distribution. 

The ratio of positive versus negative muons is a significant quantity 

which can help to discern the primary composition [7], and at high energies 
this charge ratio also reflects important features of hadronic meson production 
in cosmic ray collisions [8j. In order to obtain such muon charge information, 
we would need a way to distinguish between positive and negative high energy 
muons. Unfortunately, existing muon detectors available at shower arrays, usu- 
ally scintillators and water Cerenkov detectors |1| , are not commonly equipped 
with magnetized steel to differentiate the muon charges. Even if they were, 
the limited region of the magnetic field prevents definite determination of high 
energy muons' track curvature. 

This invites us to think of the geomagnetic field as a huge natural detector 
for muon charge information. Apparently, after being produced high in the 
atmosphere, a positively charged muon would bend east on its way down 
while a negatively charged muon would bend west, introducing an asymmetry 
into the density profile of the shower front. If their separation is large enough 
as compared with other circularly symmetric "background" deviations, it will 
be possible to distinguish the positive muons from the negative ones. 

To see such an effect, we may need a detailed simulation of inclined air show- 
ers, which keeps track of both the charges and lateral positions of numerous 
muons produced in the shower cascades. However, a simple model will be 
enough to analyze the practicability of our approach without cumbersome 
computations. We will calculate the muon lateral distribution and see its de- 
pendence on the shower zenith angle. By introducing a quantitative criterion 
to the muon energy, we can get an estimation of the condition for an unam- 
biguous geomagnetic effect. 

The article is organized as follows. We start by a general consideration of muon 
production and lateral deviations in section [2J The expressions for muon lateral 
deviations are derived explicitly in section [31 Then we present a revised Heitler 
model in section H] to complete the calculation of muon lateral distribution. 
Our criterion is applied and some typical situations are discussed. Finally in 
section \5\ we calculate the muon energy spectrum based on our model, and 
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obtain a rough image of the muon number density. We will conclude with our 
proposed approach to measure the muon charge ratio in section [61 



2 General Features 



To understand the development of an extensive air shower and the major 
processes in which muons are involved, we first introduce a Heitler model [9J, 
which describes the air shower in a simple and analytical way but captures 
the main features of both electromagnetic cascades and hadronic multiparticle 
productions. Since our interest concentrates on muon production, we will only 
review the hadronic part briefly. 

Let us assume that charged hadrons undergo multiparticle productions once 
they travel an interaction length A/, producing a new generation of N^h charged 
pions and \Nch neutral pions with equal energies, where Nch is the multiplicity 
of charged particles in the hadron-air interaction. The neutral pions decay im- 
mediately to photons, initiating electromagnetic subshowers, while the charged 
pions traverse another atmospheric depth A/ and experience multiparticle pro- 
ductions of their own. For interactions in the energy range 10 — 1000 GeV, both 
A/ and Nch can be approximated as constants [9j, where A/ ~ 120g/cm^, and 
Nch ~ 10. 

Consider a single cosmic ray nucleon with primary energy Eq incident into 
the atmosphere. After n interactions, there are N^ = {Nch)"" charged pions in 
total with equal energies 

= Eo/{^Nch 



and the shower front has traversed an atmospheric depth 

Xn = nXi. (2) 



Let us solve Eq. ([T]) for n and replace it in Eq. ([2]), so that we have a continuous 
function of X on E, 



HEq/E) 
HiNch) 



HE) = A,^^^^. (3) 



If we adopt the isothermal atmospheric density p{h) = p^e ^^^^^ with po = 
1.225 X 10~^g/cm^ and Hq = 8.4 km [T], then the atmospheric depth X can 
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be related to height H (assuming a vertical shower) by 



oo 



X{H) = / p{h)dh = Xoe 



H/ho 



(4) 



H 



where the total atmospheric depth Xq = po^o = 1030 g/cm^. Eqs. and (jl]) 
can be solved to give a relation between the pion energy E and the height H, 



We are now in a position to consider muon production and propagation in the 
air shower. Although in the above Heitler model we assumed that all charged 
pions experience hadronic interactions after they travel an atmospheric depth 
A/, they actually have certain probability of decaying to muons before they 
interact. When that happens, the muon would inherit about 1/1.27 of the pion 
energy on average [10]. Simply inserting a factor 1.27 in front of E in Eq.([5]), 
we find the relation between the production height and the muon energy. It is 
apparent that more energetic muons are produced higher in the atmosphere 
and travel a longer distance before they reach the ground. 

We should expect that the decay probability of pions in the first generations 
is extremely small, with few highest energy muons actually produced. The 
probability approaches 1 as the pion energy rapidly drops, so we set a critical 
energy ~ 10 GeV below which all pions are supposed to decay rather than 
interact [9], and the hadronic multiplication is cut off. 

To see the muon lateral distribution, we make the approximation that their 
lateral deviations from different causes can be calculated independently. Let 
us consider muons produced at a same height with equal energies. If not for 
the deviational effects, they would all have landed at a single point, i.e. the 
intersection of shower axis and the ground plane. However, due to their trans- 
verse momenta from production and multiple Coulomb scattering with air 
nuclei [11], these muons would spread isotropically in a plane perpendicular 
to the shower axis as they travel down, resulting in a circular distribution as 
a background. 

In the presence of the geomagnetic field, positive and negative muons would 
bend in opposite directions perpendicular to the shower axis. This extra geo- 
magnetic deviation splits the original landing point into two separate centers, 
one for the beam of positive muons and the other negative. Both beams ex- 
perience lateral dispersions due to their transverse momenta and multiple 
scattering, which superposes the background circular distribution onto each 
center, resulting in a distinct double-lobed feature. 




(5) 
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Conceivably, if the separation between the two centers is too small, then we can 
hardly resolve one lobe from the other. On the contrary, if we have a separation 
much larger than the radial extent of either circular distribution, hence little 
overlap between the two lobes, then we can be confident that each lobe mainly 
consists of muons with the same charge (positive or negative). In analogy to 
Rayleigh's criterion in optics, we define the condition for resolvability to be 
that the separation of the two centers exceeds twice the attenuation radius of 
each background distribution. 

In order that the separation be large, we will seek a big geomagnetic deviation, 
which implies a long muon trajectory. In this case inclined air showers are more 
favorable than vertical showers, because the total distance that a muon travels 
is magnified by a factor of sec 9, where 9 is the zenith angle of the shower axis. 
Our above analysis can be applied directly to inclined showers if we restrict 
our conclusions in the plane perpendicular to the shower axis (the transverse 
plane) instead of the ground plane. A method for projections between these 
two planes can be found in p!2]. 



3 Muon Lateral Deviations 



Let us consider muons produced at height H with energy E in an inclined 
air shower with zenith angle 9. The total length L {slant distance) of muon 
trajectory along the shower axis is therefore Hsec9. For high energy muons 
produced early in the cascade, the lateral distance from their production sites 
to the shower axis can be safely neglected, so we think of these muons as being 
produced on the shower axis and then starting to deviate from it. 

We shall calculate the muon lateral deviations in the transverse plane due to 
their transverse momenta, multiple scattering and geomagnetic bending sep- 
arately. Since we are interested in muons with sufficiently high energies, their 
energy losses during propagation are relatively small and will be neglected. 



3.1 Transverse momentum 



The charged pions from hadronic multiparticle production have a relatively 
broad distribution in their transverse momenta. When they decay to muons, 
the transverse momentum distribution is largely maintained. For convenience 
in calculation, we can approximate this distribution by a Gaussian form 
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2Pt 



(6) 



If we denote the lateral deviation due to p^ by rp, then we have ^ = where 

p ~ E'/c is the longitudinal momentum of the muon. Replacing p^ with ^ 
in Eq. ([6]), we find the radial distribution of muons in the transverse plane, 

f{r,)drv = Jjiyr^ e ^^^-P, (7) 



or in (x, y) coordinates, 



f{x,y)dxdy = - — - e ^> dxdy, 



whose standard deviation is 

.=^M)^ (9) 



Simulation of hadronic multiparticle production with QGSJet-II shows that in 
the energy range 10^ — lO^^GeV (pt) is almost a constant ~ 0.5 GeV [14j, so 

that from Eq. ©, -^^^ = ^{Pt) ~ 0.4 GeV. Adopting this value in Eq. 
gives 



400 (L/km) 
jE/GeV) 



a, = ^Z)r:.n (10) 



3.2 Multiple scattering 



Under the approximation that different lateral deviations can be considered 
independently, we will neglect the muons' initial transverse momenta when 
calculating their deviations caused by multiple scattering, assuming the muons 
to be moving along the shower axis when they are produced. 

After traversing an atmospheric depth AX, some muons will be deflected by 
an angle as a result of multiple Coulomb scattering with air nuclei. The 
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deflection angle has a distribution that is nearly Gaussian [15 



e "^^1 



'111 



The mean squared deflection angle is related to the traversed atmospheric 
depth by 



,2\ _ Es] AX 



p(3c J xo 



(12) 



where Eg = 21MeV, and the radiation length xq = 36.7 g/cm^ in air. For 
muons with energies > TeV, we have /3 ~ 1 and pc c::^ E, so that 



(0^ 



/ 0.021 Y (AX/gcm 



-2^ 



\E/Gev) 36.7 



(13) 



We can use Eq. f lT3|) and X = Xq e from isothermal atmosphere to cal- 
culate the muon lateral deviation r^. The result is also a Gaussian distribution 
hke Eq. ([8]), but with standard deviation 



(r2)V2 0.021 



72 ^/2E V 



hn 



COS 9 



cos 9 



2L + 2 



hp 
cos 9 



1/2 



.(14) 



Now the total background dispersion is the superposition of the above two 
kinds of lateral deviations. Let us "add together" these two Gaussian dis- 
tributions, which results in another Gaussian distribution with a combined 
standard deviation, 

/(x,i/) = -i^e-"^, a = + (15) 



3.3 Geomagnetic deviation 



In this calculation we leave aside the effects of both transverse momenta and 
multiple scattering. Let us decompose the geomagnetic field B into compo- 
nents -By parallel to the shower axis and B± perpendicular to it, then we define 
coordinates {x, y) in the transverse plane such that y axis is in the direction 
of B_i. Since an EAS event takes place within a small region of Earth's sur- 
face, the magnetic field B is almost a constant for a specific shower location. 
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Moreover, for very inclined air showers B± is approximately independent of 
azimuth [12]. 



Therefore, we will simply set Bj^ to a reasonable fixed value 40 /iT, and the 
muon trajectory can be well approximated by an arc in the plane containing 
the X and the shower axis. The radius of curvature is 



hence the lateral deviation 



R 



2R 



3 (^±//iT) (L/km) 
20 (E/GeV) 



m, 



(17) 



where we have expanded the bracket to first order in (^^^ , which is small 
as seen from Eq. f|T6l) . Since positive and negative muons deviate in opposite 
directions, their separation is simply twice Xg. Putting B_i = 40 fiT gives 



12 (L/km)^ 
{E/GeV) 



m. 



To see the muon lateral distribution, we now superpose the background devia- 
tion Eq. (fT5|) onto two separated centers for the opposite muon charges. Let us 
denote the muon charge ratio by i?^, so that of the total muons 

are positive and the rest ^^j^ are negative. Putting the lobe of positive muons 
on the right and the negative on the left, we finally arrive at the distribution 

1 / 1 R, {^-jf+y^\ 

f^^''^-2^ATTR/'^^TTT/^'^)' '''' 

where a = + Up is given by Eqs. & (|T^ . and s by Eq. ( ITSj) . 

Note that i?^ is found to be almost a constant around 1.268 in the energy 
range 10 — 300 GeV [16j, yet there is no systematic measurement of its value 
above that energy. Since small changes in the value of i?^ do not affect the 
validity of our approach, here we first set to 1 for symmetry and simplicity, 
so that 

We will propose an approach to measure the precise value of i?^ later. 
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Eq. (120!) is exactly the double-lob ed distribution we expected. If we take a 
to be the attenuation radius of either lobe, then our criterion of muon charge 
resolvability is s > la. Since longer muon trajectories correspond to higher 
production sites and greater muon energies, the slant distance L in Eqs. (12]), 
(|T^ and f|T8|) is a function of E implicitly. Therefore, to see the muon energy 
that meets our criterion, we will need an appropriate L-E relation to convert 
the dependence of s and o on the slant distance to that on muon energy. 



4 Revised Heitler Model 



The advantage of the Heitler model briefly discussed in section 2 is that it 
plainly demonstrates the development of extensive air showers with simplest 
calculations. However, when we consider air showers with energies several 
orders beyond TeV, some basic assumptions of the original Heitler model no 
longer hold strictly. In fact, as the hadron energy rises from a few GeV up 
to over 10^° GeV, the multiplicity of hadron-air interactions increases rapidly 
while the interaction length of the hadrons decreases by more than a half 
[T7] . Consequently, some revisions are necessary to extend the Heitler model 
to higher energies, which will result in a more realistic L-E relation than 
Eq. ©. 

Fig. [T] shows the multiplicity of charged particles Nch in pion-air interactions 
predicted by simulations using the QGSJet-II model (dashed curve). We can 
see an exponential growth of Nch with the logarithm of pion energy. This trend 
can be fitted by a function like 

iVe,(E) = Ae-^s^, (21) 



with E always expressed in GeV. A best fit gives the parameters A = 2.8 and 
a = 0.5 (sohd curve in Fig. [1]), which compares well with the simulated result. 

Similarly, we introduce another empirical formula to approximate the decrease 
of pion interaction length with energy, 

Xi{E)=B-ClgE, (22) 



where B and C are constants chosen to fit the simulated curve. Fig. [2] shows 
a comparison between our result with best parameters B = 145g/cm^, C = 
10.5 g/cm^ (solid curve) and that from QGSJet-II (dashed curve). 

To see the revised Heitler model, we only need to replace Nch and A/ in the 
original model with the empirical formulas (1211) and (I22p . More specifically. 
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Fig. 1. Multiplicity of charged particles in pion-air interactions from QGSJet-II 
(dashed curve) and the empirical formula ([2T]) with A = 2.8, a = 0.5 (solid curve). 
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Fig. 2. Pion-air interaction length from QGSJet-II (dashed curve) and the empirical 
formula ([22]) with B = 145, C = 10.5 (sohd curve). 

let us consider the nth hadronic interaction. The energy of the charged pions 



before this interaction is En-i = (|) ~ ]vT' ^^^^^ -^o is the primary energy 
of the cosmic ray particle, and A^„._i is the total number of charged pions before 
this interaction. By Eq. fl2T]) the total multiplicity of this hadronic interaction 



is M„ 



equal energies 



Ae 



alKE„ 



Er, 



n-l 



M, 



n—1 



(I)' 



producing Nn = N, 

Eq 



'n-i ■ ^Mn-i charged pions with 



(23) 



Meanwhile, the parent pion has traversed an atmospheric depth A„_i = B — 
C IgEn-i before it interacts, pushing the shower front to a greater total depth 
Xn = Xn-i + Xn-1- As usual, the multiplication of charged pions is supposed 
to continue before they reach the critical energy ^'^ ~ 10 GeV. 
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We have to use the above recursive relations repeatedly to derive an expression 
of En and X„ with regards to n only. Fortunately, this is possible due to the 
special forms of empirical formulas (^T^ and (1221) . The main steps are: 



^^^^(l-0.4343a)" ^ ^ 3 ^ ^ „ l-0.4343«r ^ ^^^^ 



2 



2 ^ l~(l-0.4343ct) 

2 

Just like what we did with the original model, we now solve Eq. (l25l) for n 
and then plug it into Eq. (126!) to obtain X as a continuous function of E. The 
result is 

^ ^ln{|A)+lgi? 



XiE) = In A) + b] + ^— (Ig E-lg Eo) , (27) 

^ ^ Va V2 ^ ^ y In 1 - 0.4343a 0.4343a ^ ^ ^ ' ^ ^ 



or, adopting the best-fit parameters, 

A-W^(-211g§^n5.4,o||^)g/c.^. ,28) 



Referring to Eq. (jl]) for the height H, and taking account of the shower zenith 
angle 6, we find the slant distance 

^ , -^0 sec 9 , , 

L = /iosecg In ° , (29) 
A(i/j 



where X(i?) is given by Eq. 

To make E the energy of the decay muon, we should insert the factor 1.27 in 
this expression, which results in exactly the L-E relation that we sought in 
the last section, 

L(E) = /^osec^ l^^l^- ^30) 



Fig. [3] compares the result from the original Heitler model (dashed curve) with 
that from the revised one (solid curve), their difference being so evident that 
our revision is justified. 

Let us set out to use this L{ET) function to complete our calculation of the 
muon lateral distribution and determine the muon energy range suitable for 
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Fig. 3. Relation between the slant distance of muon production site and the muon 
energy for air showers with Eq = 10^^ GeV and 9 = 70°, calculated with the original 
Heitler model (dashed curve) and the revised model (solid curve). 

our purpose. By replacing L in expressions ([9]), (fT4|) and (fT8|) with Eq. (!30|) . 
we render the separation s and deviation 2a in Eq. (!20|) functions of E only. 
Thus we can compare the variation of their relative magnitudes with energy. 



Fig. 4(a) shows with shower primary energy Eq = 10^^ GeV and zenith 

angle 6 = 70°. It can be seen that the geomagnetic deviation s begins to 
dominate at high energies, just as we predicted earlier. The transition point 
is s = 2a, which can be read from the graph to be at an energy of about 
10^'^ GeV. Fig. 5(a) shows the corresponding lateral distribution of high energy 



muons in the transverse plane, where positive and negative muons should form 
their own lobes respectively. As our criterion predicts, their separation begins 
to be distinguishable when muon energies are higher than 10^'^ GeV. Each 
type of muon with such energy would have little chance (< 14%) of arriving 
and being found in the other lobe, so that we can confidently distinguish the 
muon charges provided that their relative positions in the transverse plane are 
recorded with a good resolution. 

Actually 6 = 70° acts like a critical zenith angle, in which case the curves of 
2a and s nearly coincide for muon energies > TeV. For larger zenith angles, 
the separation s far exceeds 2a in the whole energy range above TeV. Typical 
examples are shown in Figs. 4(b) , 4(c) and 4(d) with zenith angles 6 = 75°, 80° 
and 85° respectively. Accordingly, the two lobes of either positive or negative 
muons can be easily recognized from each other, as seen from Fig. [51 

Note that now the lateral distribution function Eq. fl2U]) depends on the specific 
muon energy, and will be denoted by f{x, y; E). From Eq. (ITSl) . the centers of 
the two lobes become very close to each other at highest energies. Consequently 
there will be a technical upper limit to the muon energy, above which the two 
lobes are hardly distinguishable due to the limited detector resolution. Luckily, 
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Fig. 4. Dependence of the separation s (solid curve) and attenuation radius 2a 
(dashed curve) in Eq. (j20p on the muon energy, assuming air showers with pri- 
mary energy = 10^^ GeV and at different zenith angles 61 = 70°, 75°, 80°, 85° 
respectively. 



as Fig. [8] will show, the number of muons produced at such highest energies is 
virtually zero, hence this energy range will be excluded from our consideration 
of muon production and deviation. 



Note also that with the growth of zenith angle, the separation between positive 
and negative muons becomes larger, so that the demand for detector resolu- 
tion is more accessible. Therefore, in order to detect an unambiguous geomag- 
netic separation between the positive and the negative high energy muons, 
we should focus our attention on EAS events with zenith angles > 70°. On 
the other hand, for almost horizontal air showers with a nonzero azimuth, the 
two lobes of muons would travel very different distances to reach the ground, 
resulting in extra asymmetry in the positive and negative muon distributions 
|12j . Therefore, we suggest an optimum zenith range 75° < 6 < 85° for our 
approach, which contributes over 17% of the half total solid angle. 
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(c) e = 80° 



(d) e = 85° 



Fig. 5. Lateral distribution of high energy muons in the transverse plane as a 
function of x coordinate with y = 0, for fixed muon energy around 10^-^ GeV m 
air showers with primary energy Eq = 10^^ GeV and at different zenith angles 
(9 = 70°, 75°, 80°, 85° respectively. 

5 Muon Number Density 



We note that the lateral distribution we calculated in the early section is only 
for muons with fixed energies. To see the actual muon number density in an 
extensive air shower, we need to know the muon energy spectrum, of which 
we can get a rough estimate through the same Heitler model. 

We start by considering the number of charged pion produced in the develop- 
ment of an air shower. From Eq. ( l23i) we have 



or, eliminating n by Eq. ( l25l) and adopting best-fit values of the parameters, 

^ ^^/ 2.870 + lgEA"-° (32) 
^ ' E l 2.870 + lgEoy ^ ' 



In this expression we have assumed the proliferation of charged pions to be 
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Fig. 6. Estimated probability of pion decay at high energies, calculated for air show- 
ers with primary energy Eq = 10^^ GeV and zenith angle 9 = 80°. 

a continuous process, during which the number of charged pions grows while 
their energy drops simultaneously. We can calculate the number of pions pro- 
duced in the energy range E to E — dE by taking the derivative of N{E), 




As mentioned in section 2, the parent pions of these newly produced ones 
actually have a small chance of decaying to muons instead of interacting. We 
can estimate this probability by taking the ratio of the interaction length to 
the decay mean free path, 

P(E) ^ ^^(-^) = ^0"^^^ ME) ^ 115 145- 10.5 IgE 

^ '~ c-iT^±p{H) CT^±E X{E) cos 9 E X{E)cos9 ' ^ ' 



where we have used the relation between the atmospheric depth X{H) and the 
isothermal atmospheric density p{H). Since the probability of pion decay is 
very small with energies above 10^ GeV (see Fig. [6]), we can safely neglect the 
charged pions lost by decaying in the Heitler model before the energy reaches 
too low. Thereby the number of high energy pions that decay in the energy 
range E to E — dE is 

N^^(E),E^P(E)^^[yjp),E, (35) 



where the multiplicity of charged particles Nch{E) is given by Eq. (l2Ti) . 

The muons from these pion decays typically lose several tens of GeV's energy 
to ionization before reaching the ground, which is negligible compared to their 
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Fig. 7. The energy spectrum of high energy muons in air showers with primary 
energy Eq = 10^^ GeV and zenith angle 6 = 80°, expected by the revised Heitler 
model. 

kinetic energies well above TeV [11]. Hence their energy spectrum D^{E) 
remains virtually the same before they reach the ground. Allowing for an 
energy fraction factor in pion decay, we have 

D^{E) = 1.27 Ndec{l-27E). (36) 

Fig. [7] shows the calculated result for D^{E) with shower primary energy 
Eq = 10^^ GeV and zenith angle 6 = 80°, exhibiting a perfect power-law 
oc E"' with index 7 ^ —3. 

Let us integrate Eq. (l36l) to see the number of muons with energies > E, 

Eo 

N^{E) = J D^{E)dE, (37) 

E 

where Eq is the energy of the primary cosmic ray particle. Moreover, with 
the muon energy spectrum, we can calculate the "integral" number density 
of high energy muons in the transverse plane. Referring to Eq. (!20!) for the 
lateral distribution of muons with energy E, let us multiply it by the spectrum 
D^{E) and integrate. 

El 

g{x,y) = J f{x,y;E)D^{E)dE, (38) 

Eb 

where we have chosen the lower bound of integration to be E^ = 10^ GeV, and 
the upper bound Ei ?a 10^'^ GeV is given by N^{Ei) = 1, above which energy 
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Fig. 8. Expected number of high energy muons in the revised Heitler model for air 
showers with Eq = 10^^ GeV and 6* = 80°. 
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Fig. 9. Contour plot: number density of high energy muons in the transverse plane 
for air showers with primary energy 

^0 = 10^^ GeV and zenith angle 9 = 80°, 
calculated with the revised Heitler model. 

there is no muon produced. The numerically calculated muon number density 
g{x, y) for 6 = 80° is shown in Fig. [9l which can be compared with simulation 
or experiment results. 

As expected, the number of high energy muons is very limited (Fig. [8]). So 
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the detector array will actually record several high energy muons arriving 
at discrete lateral positions, instead of a continuous density distribution in 
the transverse plane. Nevertheless, once we have determined the shower core 
position, our calculated double-lobed distribution (Fig. ^ predicts that the 
high energy muons arriving at each side of the shower core most probably 
(> 99%) belong to the lobe on the same side, carrying the corresponding 
charge. The shower core can be located by analyzing the density distribution 
and energy deposits of plenty lower energy secondary particles ^8j. Since we 
know which side for which type of muon charge from the direction of local 
geomagnetic field, in this way we can confidently identify these high energy 
muons with their own charges. 

Our final purpose is to obtain the charge ratio for such high energy muons, 
which will give us valuable information on the composition of primary cosmic 
ray particle and the understanding of hadro-production processes. However, 
we need to point out here that it is impractical to measure the charge ratio 
on an event-by-event basis. Because of the stochastic nature of extensive air 
showers, the fluctuation in muon numbers will be large compared to the limited 
total number of high energy muons in a single event. Also, not all the high 
energy muons in a particular event will be detected, since the array detectors 
do not cover the whole region. Therefore, both positive and negative muon 
numbers should be collected from a large number of inclined EAS events, so 
that we can get a statistically sound value of for high energy muons. 



6 Summary 

In this article we analyzed the possibility of obtaining the charge information 
of high energy muons in very inclined extensive air showers. We have demon- 
strated that positive and negative high energy muons in sufficiently inclined 
air showers can be distinguished from each other through their opposite geo- 
magnetic deviations in the transverse plane. We developed a revised Heitler 
model to calculate this distinct double-lobed distribution, and studied the con- 
dition for the two lobes of either positive or negative muons to be separable 
with confidence. From our criterion of resolvability, we concluded that a zenith 
angle 75° < 9 < 85° will be most suitable for our approach. 

There are already some results from full air shower simulations that take into 
account the geomagnetic effect on muon propagation [12] [IS]. They illus- 
trated remarkable double-lobed muon lateral density profile in very inclined 
air showers, which is in agreement with our expectation qualitatively. However, 
no present study has fully considered the high energy part of muon content, 
which can be used to compare with our results. Thus we would hke to pro- 
pose future simulations of very inclined extensive air showers that focus on 
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the behavior of high energy muons. They also have to keep track of the muon 
charges and the relation to their lateral positions. 

For our method to be applicable in measuring the high energy muon charges, 
there are some requirements on the shower detector performance. First, the 
muon detectors should be able to measure muon energies up to over 10^ — 
10^ GeV, so that we can sort out those muons with energies high enough 
for our purpose. Moreover, since the high energy muons arrive very near the 
shower core, we would have to look for EAS events whose shower cores are 
inside the coverage of the detector array. Besides, to distinguish the muons' 
charges from their lateral positions, a detector array is expected to have a high 
resolution in the shower transverse plane. 

As far as we are concerned, there are still some technical limitations to the 
above requirements. For example, the detectors near the shower core are usu- 
ally saturated with signals from lower energy muons [1]. Therefore, we suggest 
to employ muon detectors that are especially aimed at detecting high energy 
muons in the shower array. They should have a high threshold energy up to 
10'* GeV, so as to avoid unwanted signals from low energy muons. Or these 
detectors may be extensions of existing ones like water Cerenkov detectors, 
but with special triggers to sort out high energy signals. 

Our method not only works for air showers with primary energies over the "an- 
kle" (~ 10^° GeV), it can also be applied to study cosmic rays in the "knee" 
region. In such cases, the primary energy of the cosmic ray particle is around 
10^ GeV, several orders lower than the examples we studied above. Accord- 
ingly, the energies of muons produced by early generation pions range from 
10^ GeV to 10^ GeV. Produced higher in the atmosphere with lower energies, 
these muons experience larger lateral deviations, and are more susceptible to 
the geomagnetic bending. Thus we shall have an even better condition to use 
our method to obtain the muon charge information. What is more, for such 
energies below 1000 GeV, there are alternative ways to distinguish the muon 
charges [20||8] . By comparing the results, these experiments can serve as a test 
to validate or rule out the possible application of our approach. 
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